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Abstract. We investigate the cylindrical class of Newman Tamburino equations in 
the presence of an aligned Maxwell field. After obtaining a complete integration of the 
field equations we look at the possible vacuum limits and we examine the symmetries 
of the general solution. 



1. Introduction 

In a previous work we examined the generalization of the Newman Tamburino metrics 
in the presence of an ahgned Maxwell field [1]. It was shown there that the so called 
'spherical' class did not admit any solutions and that consistency of the field equations 
therefore required the 'cylindrical' condition \p\ = \a\. The existence of solutions in the 
cylindrical class was left however as an open question. 

In the present paper we obtain a complete integration of the Einstein Maxwell 
field equations for this problem. We first show that the Newman-Penrose, Bianchi and 
Maxwell equations form an integrable system and then proceed to integrate the first 
Cartan structure equations. After obtaining the general solution we look at the limiting 
case, where the charge Q of the Maxwell field goes to zero. Surprisingly in this limit, 
we do not recover the empty space metric obtained by Newman and Tamburino, but 
rather the special case admitting two Killing vectors. 

In §2 we describe the choice of tetrad and the integration method. The resulting 
metric is presented at the end of this section. In §3 we discuss the vacuum limit of the 
metric and its symmetry properties. 



2. Cylindrical class 

We follow the notations and conventions of our previous paper [1]. The cylindrical 
class of Newman Tamburino metrics is characterized by the existence of a principal null 
direction k of the Weyl tensor which is hypersurface orthogonal and geodesic, 

^0 = 0, (1) 



K = 0, 



(2) 
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p-p = 0, (3) 

but has nonvanishing shear and divergence, with the spin coefficients p and a related 
by 

- W\' = 0. (4) 

We assume now that a non-null Maxwell field is present, which is aligned in the sense 
that k is also a principal null vector of the Maxwell tensor : 

$0 = 0. (5) 

As mentioned in |1] we can, by means of a rotation, a boost and a null rotation, set up 
the complex null tetrad (m, m, 1, k) such that it is parallelly propagated along k : 

K = e = 7i = (6) 

and such that, additionally, 

T = a + p. (7) 

The NP equations guarantee then that D{a/a) = and hence, by (jlj), a spatial 
rotation exists such that 

a = p. (8) 

The remaining tetrad freedom consists now of boosts and null rotations with parameters 
A and B respectively satisfying : 

DA = 6A = 6A = = DB. (9) 

This allows us to rewrite the NP, Maxwell and Bianchi equations determining the 
evolution along the geodesic rays as follows : 

Dp=2p^ 

D a = {a + P) p, 

D X = {X + p) p, 

Dp = {X + p)p + ^2 + ^R, (10) 

D u = {a + p) p + {a + p) X + + ^2'^, 
DP = {a + P)p + ^u 

D'j = aa + 2ap + pp + '^i- ^R+^i'^, 



while 



and 



D<l>i = 2p$i, (11) 
L)^i = 4p^i, (12) 

= 2(a + /?)$i -p$2, (13) 
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(5^1 = 2(2 « + 3/3)^1 - (3^2 + 2 $1 $7) P- (14) 

Next we calculate [6, D] from which we eliminate D \E'2 by (2 i?9 — 3 B^) : 

_ 3 

-^i5p + p6-^i + 2p^<l>i<l>i - (a + /5 + 3a + p^i - -^1^ = 0. 

The last equation, together with [S,D] (p^&i), gives us expressions for 6 p and : 



(15) 



5p =\(5^^-2^^p' + ^{9{a + P)-5{a + P))p + ^^, 

+ -p^2-2p$i$7. (16) 
From (2 i?9 — 3 B^) we find the evolution of \I^2 along the geodesic rays : 

We can now simplify NPu and its complex conjugate by means of [6, D] p, and eliminate 
\I^2 from both. This results in 

2 p 

and 

2(a-/3 + a-/3)p + ^i+^=0. (18) 

From [6, D] $1 and the second Maxwell's equation we get expressions for D $2 and 
which we can use to solve [6, D] $1 for $2 : 

1 / - 2^1 -¥r\ 

D$2 = - \ 3a + f3-a + f3 + -j $i + p$2, 

=l(^3« + /3-a + ^ + ^^^i— ^^$1, (19) 

$2 =-{a + 3f3 + a-(3- — ]—. 

2 V P J P 

In the next step we use (fT^ [T5l [T6|) for evaluating [5, and [6, 6]p. Herewith 
we obtain 

S{p + p) = Sa + 26a-6a + (a - P + 3 p - 4 ^^^^^ ^ (3 

-(2« + 3«-4/3 + 2^)/? + 2(p-p + 2^^)p (20) 



— 2a — aa + a — [a — a) — 

P 
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and 

5{p + p) = Sa-5a + 25a + + 2a - 2/3 - 4^^^^ + 3y ) p 



+ (4a-a-3(3-(3 + ^yp-2^p-p-2^^y (21) 



2 - ^1 

— 3a; — 3aa — 2a — 
P 

which allows us to find the 6 derivative of ffTSl) as 



5a-5a = (/i-/i)p+^(3a^ + /3^ + a^-9/9^)+aa-2/5/3 

- 3-1 1 / - ^1 

+ a/3 + -a/3--a/5+- a + 3/3 + a + ll/3 — (22) 

2 2 8 V P J P ^ ^ 

+ ^(^a-P + 9a-P + 3^^ ^ + (^2 (/3 - a) ^ - 1^ <l>i$7. 

Now follows a crucial step, in the form of two beautiful factorizations. Expressing 
that the right hand sides of (!20!) — (pT!) and (1221) + (!22|) are and eliminating from both 
of these a by means of equation (fT8|) . we find that: 

(2/3p-2pa-^i) (p2$i$r*7 + 3p2$i$T^i +^^1^) =0. (23) 

and 

("^-^i) (2/3p-2pa-^i) (2p2$i$7+^i^) =0. (24) 

It can easily be seen from the above results that the common factor 2/3p — 2pa — \&i 
has to be zero. If not, then by would have to be real, after which the second 

factor of ([23D would factorize as times a positive definite expression. We therefore 
conclude that 

/5 = a + - — . 
2 p 

Next we show that the Ricci scalar (or, equivalently, the cosmological term) is zero: 
R = 0. In order to do so, we first rewrite NP12 as 

/ 1 _ 1 ^1 -W[\ 
(5a = 5a + (A — u)p + 2(a — a)a+ 2a aH — 

_ V 2 4 p ; p 

rt, 

2 p 24 ' 

after which the fourth Maxwell's equation reads 

+2 (a + (5a - aa) — -2p$i-A$i = 0. (25) 
P 
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If we calculate the propagation of this expression along the geodesic rays, eliminate 
the terms D A $i and D6a hj means of the commutators [A,D]$i and [5, and 
simplify the result by equation ( l25i) and NPiq, we get 

- i 11a + a + 2 — ] — - - ( a + Ua + 2 — ] — + 6 (a^ - a^) 
2\ p J p 2\ p J p ^ 

-2{p- p + 2-f -2^) p + 2(6a-Sa) - Ir = 0. (26) 

6 

Adding (126!) and its complex conjugate results in i? = 0. 

From [6, D] a, NPiq and {3 Bg + 2 Biq) we now see that respectively 

/ 13\E'i+l^"\ 

D6a = Ap6a + (X - p) p^ + { a + a -\ + a^i - p <l>i$i, 

V 8 p J 

1 / ^\ ^1 1 , -X ^ / 

An = — \ a + 5a -\ (3a — 5 a) 1-(37 — 7 — 2u)p 

4V P J P ^ P 



+ - Aaa - a^ + 25 a - $i$i 



and 



A^i = ^ (3 ^1 - 5 % + ^ ((3 a - 25 a) ^1 + (3 a + 7 a) ^1) ^ 
8 p'^ 4 

+ i25a-Sa^-aa + 2a'^- <l>i$7) — + 2 (2 7 - 7 - p) ^1 

P 

+ 2p^3 - 2 (a + a) <l>i$T- 
Then we rewrite the Bianchi identity (3 -B7 + 2 Biq) as an equation for D \l/3 : 



D^, = 2^+[^a-a--^]^ 

p3 V ^ P J P^ 

/ „ 3al^"\ 

+ + 2 5a - 4 a a - 2 pp - $i$i h 2p^3. 

V P J P 



(27) 



We now try to find an expression for 5 a. In order to do so, we evaluate the 5 
derivative of equation fl26|) . Eliminating from the result 5 7 and 57 by A^Pis and the 
complex conjugate of NPis, we get 



4 V J p-^ p-^ 

+2 (3 (2^1 - ^) a + (2^ - a) + (3 a - 11 a) 

+ (4 5a + 14a^-2aa - — 

P _ 

- (4 (5 a + 4 - 14a a + 22 - $1 $7) — 

+ (7 - 7 - 4 p) ^1 + (3 (7 - 7) + 4 p) ^ + 4 (z/ - z/) p2 

+4 ( A a - A a + (7 - 7 - 2 p) a + (3 (7 - 7) + 2 p) a + - ^3) P 
+8 [{a- a)5a- aa^ - a^ + 2 a a^) + 4 (a - a) = 0. (2^ 



Complete integration of the aligned Newman Tamhurino Maxwell solutions 6 

Herewith the [5, A]p-relation yields an expression for while from the [5, S\ a-relation, 
we get an expression for 55 a. This then enables us to simplify [5, A] $i, which results in 
an expression for A a (note that the expressions we obtain all contain terms with factors 
5a and/or 5 a va. the right hand side). Substituting the latter into (128|1 . we eventually 
obtain : 

p2 5a = (7^?+ (8 (6a- a) p- 11^)^?) 

+ — f 16 p - 32 (2 a - a) a p2 - 4 (5 a - 9 a) - \ 
16 V / 

- (W^ - ^3) p^ + p^p^ + (-3 a^Wi + 2a aWi - a^Wl) 

(12 aWi^ -SQaWf^ p-3Wf. 

Now we can write A^Pi5 and NPi^ as expressions for 5 7 and 5'y, after which we can 
solve (1261) for 7 : 

^ = ^ + + 2 + s 3 • 

p p^ 8 p'^ 

Eliminating (5\I'3 from [5, A] $1, we can also solve [5, A] \l'i for \E'4 : 

-'-^%-3(2a^ 

8 p ; p5 1^ 2 p 





8 p3 64 p4 ) p2 

2 p ; p • 

As in paper [2] , we proceed by isolating the p-dependence of all variables. Equations 
( ITT]) and ( fT2l) allow us to introduce variables 0i and V'l, respectively defined by: 

$1 = (pip, D(j)i = 0, 

^i=^ip2, D^i = 0. 

We also define a variable L as 
L = log IpI . 

Since 

D{a/p) = ^ij,DL, 

we see that 

a = ^tpiLp + ap, Da = 0. (29) 

This enables us to put a equal to zero by a null rotation with parameter B = 
{jipi L p — a) / p (as this is compatible with the condition D B = Qj). 
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Comparing the derivatives of a as in fl29p with those we aheady obtained before, 
we can extract expressions for A, yU and v : 



\ = -[L^ + ^L + '^\p^'i-^{?>L + 1) p^i^i - ^ pL0i0i 



4 

1 fl 



^3 



) + ( ^V^i^ + 0101 



+ i7V--i(r^ + i)M^,^i-i((i + |)L-2)^. 



Herewith we can integrate the D \l/3-equation ( 1271) . which gives 

*3 = 



^ (L + 3) L^IjI - ^L^iVi - ^^^10101 + ^3 j p\ D^Ijs = 0. 

We now come to a more subtle part of the integration. First we introduce a help 
variable Q = 6 p/ D p and define two new operators ei and 62 having the property that 
ci p = 62 p = and Deix = De2X = for all x obeying D x = 0. One easily sees from 
the [6, -commutator that this can be achieved by putting 6 = ei+pe2+^D and 
6 = pe2 — ei + QD. The new commutators [ei, D] and [e2, D] are then given by 

[e^,D] = {[6,D]-[lD])/2, 

[e2,D] = {[6,D] + [lD])/{2p). 

The expressions for the commutators [ei, A] and [e2, A] can be derived from [S,A] and 
[6, A] . They read 

, 2 



[ei,A] 
[e2, A] 



Y^pL [tpi + i^iY ei, 



1/2 2\ 1 

^ (^^1 - V"! j ei + 2 7 62 - - 



+ -L + 2 ) 



(L + 3) ^1 - ^V^i j L^i+2(L + 1)0101 



pe2 



1 

^ 64 



(L + 7) - (2L2 + 5L - 3) ^I'V'i 
(^(^2L^ + + + 4 (2L2 + L - 1) 0107^ 



{L" + 2L-2) V^i0i0i + (L + 1) V^i' - 32 ( ^ - 1 ) V^; 



^1 
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Finally, from [S, 5] we get an expression for [e2, ei] : 

1 

"i 



1 

[e2,ei] = -- (^/^i + -01) ei. 



Note also that 

ei = — ei and ^2 — ^2- (30) 

At this stage we have expressions for all derivatives of 7, p (thus of L), 0i, ■^i and 
■03 ( can be found from NPiq) . 

Next we integrate the D 7-equation and obtain 

7 = ^0 + ^ (^1 +-01) pL"^^! + ^pLipi^ + Q^i^ + ^4>i4>i^ P, 

in which we can put go equal to zero by a boost. To demonstrate this we first look at 
the derivatives of 0i, which are all zero, except for A0i, which equals — 4 0i 0i gQ. As a 
boost transforms 0i into (pi/ A, we can put go to zero by choosing AlogA = —2gQ. This 
also allows us to write 0i as Q f, with Q a non zero constant, and with / on the unit 
circle (/ = !//). 

Now we write ipi as ^1 = {U + V)/2 with U real and V imaginary. If we then look 
at the derivatives of /, 0i, U and V we see that 

df = e^f + ^f=l^{A + LpU), 

d0i = 62 01 + A 01 = {A + LpU), 

dt/ = e2t/ + A^ = ^^^^ i^__^(4 + Lpt/), 

dy = e2y + Al^ = ^(4 + Lp[/), (31) 

showing that /, 0i, U and \^ are all functionally dependent. Notice that U cannot be 
constant, as then V would be real : 

1/2 TJ2 

e2C/=^-^-4e (32) 

This permits us to use f/ as a coordinate, but we prefer to write U = U [x] (and hence 
also V = V{x), f = f{x) and 0i = 0i(a;)), with x to be specified later. 

At this stage the only further information comes from introducing appropriate 
coordinates. First we construct the basis one forms Q}, Q?', fi^, 11^, dual to Ci, 62, A, 
D. 

The Cartan structure equations imply dQ^ = and thus — <\u. 
Since dp — ^ pQ.^ -\- 2 p^ Vt^ we have 



4 dp A p dw 
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Next notice that dU = Vi^ + VL^ and = ^VVL^ + 9? , where, by ^ 
and fl32|) . 62?/ and cannot be 0. It follows that Vt^ can be written as a linear 
combination of df/ and du. As AU = L pU {2 V"^ — — 16 Q^) /16, we can write 

= S{x) dx - ^^d«, 

with S{x) to be specified below. There remains which is general, = B du + C dp + 
H dx + J dy, but where we can choose C = by a transformation of the y-coordinate. 
From fl30l) we see that 

J = -J, 
S —— s ^ 
H = -H. 

The tetrad basis vectors are then given by 

D = 2p'^^. (36) 

Acting now with the commutators on the coordinates u, p, x or y, we get 
DJ = DS = DH = = ei5, 
AS = ^LpUe2S, 

DB = ^pV. 

It can then easily be seen that 

B = ^LV + Bo, DBo = 0. 

From D Bq = = DH = DJ one deduces the existence of a new y-coordinate, as a 
function of u, x and y, such that H becomes 0. 

From [e2, ei]y we find the following expression for 62 </ : 

^ Jf/ dlogJ 5f/ 
62 J = or — ; = , 

4 da; 4 ' 

from which we see that J = jQ{x)Ji{u,y). Defining a new y coordinate as a function 
of u and y and absorbing the dw-part in a new S-coefficient, one can assume Ji = 1. 
Herewith J = J{x) and we have that ei J = and A J = L pU'^J/16. 

From [ei,A]?/ we find that ciBq = 0, which means that we can write Bq = 
J Bi{x,u), after which [e2. A]?/ results in 

62^1 = -f/y/(4J). (37) 
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Herewith Bi can be decomposed as Bi = B2{x) + B^{u). A final w-dependent y- 
translation allows to transform VLi into Vti = {VL/A + JB2{x)) du + Jdy, i.e. one can 
assume Bi = Bi{x). 

We now fix S{x) such that we can integrate (!37j) : As 62 V = UV/A and 62 J = UJ/i, 
the choice 

S{x)=4/{xU) (38) 

leads to 

V 

Bi = -—hgx. 

From the 62 l^-equation, we see that V = ax, with Qa constant {Q is assumed 
to be non zero). We still need a solution for U. From the expressions for DU, AU, 
ei U and 62 U it is easy to see that U = U{x) with f/ (x) determined by 

(9f/ _ _16Q2 + 32Q2^V + ?72 
which integrates to 



(39) 



From ei, 62, A and D applied to (^ip3 + /2 we see that the partial derivatives of ip 
(the real part of ip^ divided by U) are given by : 

dtp ^ dip dip 
dy du dp ' 

dip _U'^ + 16 (1 - 14xV)g2 _ 64 V 
dx 32 X 

Using the equation (HOl) we can solve these differential equations for ip : 
1 

ip = - [a c?logx + 8 C2 — 15 a^x"^) — , with ci, C2 constant. 
8 ^ ^ x^ 

The dual basis now becomes 

/I \ 2 1 

fi^ = I - (ilQax — U) L — 2 IQaxlogx ] du H — t^— dx H — IxJody, 
\8 J xUp 2 

fi^ = ( -- UlQax + U)L + 2 IQaxlogx ] du + dx - -IxJody, 
\ 8 / xUp 2 

= du, 

~ 1 



= 7-^dp - ^ dx - ^—^JoQx'^ady + 
2 p'^ 2x p 2 



^(2 + L)LU'-2iP+ 
64 



(^2 (logx - 2) x^a^ + ^ (1 + (41og3; - 1) x^a^) L - ^x^a^L^^ 



du, 



which shows that we can absorb Jq in y. 
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The metric becomes 



1 



px 



1 



(2L — 41ogx + 1) Qx^ady -dp H dx ] du -\ — x^dy"^ 



+ 



2L^a^x^ - + - 2)a^x^^ L + 2 (^^ - 2a^x''^ hgx (41) 



4 C2 a^x^ 



+ 



2 {cl - x2 - a^x^) p^Q^ 



dx^ 



with the Maxwell field given by 
$0 = 0, 

$1 = Q/P, 



$2 = 1 {LU + 4(L + 1) jgax) g/p, 



(42) 
(43) 

(44) 



where 



/ 



2 av/c? - (a2x2 + l)x2 + (2 a^x^ + 1)/ 



V4a2c? + 1 
is on the unit circle, as required. 

A more elegant expression for the metric dHJ, preserving the obvious symmetries 
of the solution, is obtained by a coordinate transformation p ^ rx"^, which results in 



+ 



2{c\ - x2 - a2x4)r2x4Q2 

a2(2L + l)2g23;2 (Lc2-4C2)g^ 

2 a;2 



dx^ + a(2L + l)QxM?/ 



dw^ + —dy", 



dr I du 



in which L = logr = logp — 2 logx. 



3. Vacuum limit 



Neither u nor y appear in the components of the Riemann tensor \Efi, '^2^ ^3, ^4 and 

$0, $1, '^'2 (BZD-dllD: 

^ _ ([/ + 4jgax) p2 
^ ~ 2 ' 

v^2 = + ((L + l)t/ + 4(L + 2)/gax) , 

o 



(f/ + 4/gax) p2 
128f7 



(3L + 4)t/^L + 24(L2 + 3L + l)IQxaU^ 



-16 ( 2L + (3L^ + 14L + 9)x^a^ - 11^^41^') g^f; 



x^ 



^4 = + 4/Q«a:) p2 ^ ^ ^ r.g^2 ^ ^ l2)IQxaU^ 

256(7 L 
+ (L + l)t/3L - 16 (21 + (L^ + 7L + 4)x2a2 - 41ogxcf + 8c2 \ ^2^" 
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and one can easily see that also in the canonical Petrov type I tetrad precisely two 
functionally independent functions will remain in the curvature components. Therefore 
the metric will admit exactly two Killing vectors [3]. The vacuum limit of our metric 
then can not be the general one obtained by Newman and Tamburino, as the latter has 
only one Killing vector. In fact the {Q = 0)-limit of (HTl) is given by 



which, after the coordinate transformations x Cox/-\/2, p — > — l/(2r), ?/ — 2y/co, 
equals (26.23) of ^. 

4. Conclusion 

We have obtained the general solution of the 'aligned Newman- Tamburino-Maxwell' 
problem. These are the algebraically general metrics possessing hypersurface orthogonal 
geodesic rays, with non vanishing shear and divergence, in the presence of an aligned 
Maxwell field. We have shown that the solution necessarily satisfies the cylindrical 
condition \p\ = \a\ and has exactly two Killing vectors, which implies that the vacuum 
limit is only a subfamily of the general cylindrical Newman Tamburino metrics. 
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